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is discussed. 

PACS numbers: PACS number(s): 03.65.Bz, 42.50.Dv 



I. INTRODUCTION 

A number of different phase space distribution func- 
tions has been introduced and investigated over the years, 
the Wigner distribution function 0] being the most fa- 
mous and the most known of all. The phase space dis- 
tributions from the Glauber-Cahill s-parameterized class 
of quasi-distributions j^] that contains the Wigner func- 
tion, the Glauber-Sudarshan P-representation and 
the Husimi or the Q-representation Q , have been widely 
used as useful and powerful phase space tools. 

It is the purpose of this paper to show that a differ- 
ent class of phase space distribution functions, connected 
with a lesser-known Kirkwood distribution function 
can be useful in the investigation of quantum interfer- 
ence in phase space. The Kirkwood distribution was 
proposed just a year after Wigner introduced his func- 
tion and like the Wigner distribution was firstly used in 
quantum statistics and thermodynamics. The Kirkwood 
distribution has been rediscovered by Rihaczek |^ for 
use in the theory of time-frequency analysis of classical 
signals. Operator bases of Kirkwood and Wigner distri- 
bution functions have been examined by Englert |^ . Zak 
1^ studied the Rihaczek function in the quasimomentum- 
quasicoordinate representation. 

In this paper we shall study a new class of phase space 
distribution functions connected with the Kirkwood- 
Rihaczek distribution function. The main advantage of 
the functions from this class is that they all lead to the 
correct position and momentum marginals. We shall 
show how to construct a coherent state representation of 
such functions and which operator ordering corresponds 
to the Kirkwood-Rihaczek distribution and their gener- 
alizations. A connection of such phase space functions 
with squeezed states will be established. 

In Section II we present the basic definitions, prop- 
erties, and differences of the Wigner and Kirkwood dis- 
tributions. Section [II is devoted to the problem of op- 



that a full quantum state reconstruction with the help 
of the Kirkwood-Rihaczek function is possible. Sections 
^ and VI are devoted to the quantum interference in 
phase space. Using simple examples, we present various 
properties of the Kirkwood-Rihaczek functions and we 
discuss similarities and differences of such phase space 
quasi-distributions with the Wigner function. 



II. PHASE SPACE AND QUANTUM 
MARGINALS 



A. The Wigner distribution function 

In 1932 Wigner introduced a phase-space distribution 
function 



erator ordering and its relation to various phase space 
distributions. In Section IV the definition of general- 
ized Kirkwood-Rihaczek distribution function is intro- 
duced and it is shown that Kirkwood-Rihaczek func- 
tion fully characterize the quantum state. This means 



(1) 



which fulfills the fundamental requirement for a joint 
probability distribution in phase space, i.e. when inte- 
grated over q or p, it gives marginal probabilities: |5'((z)P 
and \^{p)?, where = J dqe"^ '^{q) is the 

Fourier transform of '!'((?). 

As shown by Wigner, the function given by Eq. (|l|) is 
in general not positive. However, under simple and rea- 
sonable physical assumptions it is unique. Thus, when 
one demands that: (i) a phase space distribution P{q,p) 
is real, (ii) bilinear in 5", (iii) gives the correct marginals, 
and (iv) its dynamical evolution reproduces the Liou- 
ville equation in the classical limit, then the distribution 
P{q,p) is unique and is the Wigner function | p^ . 

After Wigner's work many different phase space distri- 
bution functions have been introduced and investigated. 
There is a very rich literature devoted to applications of 
the Wigner function and other various quasi-distributions 
in quantum optics (see e.g. [HI). 

The simplest example of a phase space distribution 
which does not satisfy the four assumptions leading to 
the Wigner function, but reproduces the right marginal 
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in an explicit way, is the distribution function: 



(2) 



Clearly, this distribution is not bilinear in 5* and as a 
result does not satisfy the requirements of the Wigner's 
uniqueness theorem. Another distribution function sim- 
ilar to that given by Eq. (|^) can be guessed taking for- 
mally a square-root of this expression. As a result, up to 
an arbitrary phase, we have 



P{q,p)^\^{q)\ \^{p)\. 



(3) 



Note that this distribution function contains no infor- 
mation about the phase of the wave function. A simple 
insertion of an additional phase factor Lp{q, p) to the wave 
functions leads to the expression: 



P{q,p) - e^'^^^'P) ^*ip) 



(4) 



that defines a class of bilinear but complex distribution 
functions. An example of such a distribution has been 
proposed just one year after Wigner introduced his fa- 
mous distribution function. 



B. The Kirkwood-Rihaczek distribution function 

In 1933 Kirkwood introduced a phase space distribu- 
tion which, according to his description, "... differs but 
little from the Wigner function" . The Kirkwood func- 
tion is defined as follows: 

K{q,p) = ^/ de*(9)e'^**(0 = 



(5) 



In 1968, this function was rediscovered by Rihaczek in 
the context of a signal energy distribution in time and 
frequency (for a review of time-frequency distributions 
see nil). 

The real part of the Kirkwood-Rihaczek distribution 



KRc{q,p) = Re 'i'iq)e 



(6) 



is closely related to a quantum mechanical phase space 
distribution introduced by Margenau and Hill p^ . 

We clearly see that the Kirkwood-Rihaczek (K-R) dis- 
tribution has the correct marginal properties: 



K{q,p)dp = \^iq)\' 



K{q,p)dq = — 



and it is normalized. 



K{q,p) dqdp — 1, 



(7) 



(8) 



as a consequence of the normalization of the wave func- 
tion Note that the contribution of the imaginary 
part of the K-R distribution to the marginals is identi- 
cally zero. Hence, we shall often focus on the real part. 

Another condition easy to find is that the absolute 
square of K{q,p) has the form similar to the Eq.(||), i.e. 

\K(.q,p)\'=j^\^iq)\'\^ipr, (9) 

which indicates that K(q,p) is a square-integrable func- 
tion, with 



K{q,p) P dqdp = 



1 



{2T:h) ' 



(10) 



The absolute square of the K-R function, Eq. (^, has 
a simple physical interpretation. It is just proportional 
to the product of the probabilities in configuration and 
momentum representations. 

The dynamical free evolution of a particle with mass m 
of the K-R distribution function is given by the following 
equation 

dtK{q,p,t) + Ld^K{q,p,t) = ^dlK{q,p,t) , (11) 
which can be also written in a form: 



Kiq,p,t) 



^>Kiq,p,0). 



We see from this formula, that the free evolution of the 
K-R function is a superposition of the free Schrodinger 
diffusion and of a classical boost to a moving frame. The 
K-R distribution function is bilinear and has the cor- 
rect marginal properties but, in contrast with the Wigner 
function, it is not real nor does its free evolution satisfy 
the classical Liouville equation. As we shall see in de- 
tail later, it is also not well-behaved under rotations of 
the {p, q) coordinate system, and this implies that such a 
function can not be measured by tomography methods. 

The simplicity of the definition, Eq. (|g), indicates that 
it is relatively easy to evaluate the K-R distribution func- 
tion even for systems for which an analytical formula of 
the Wigner function is not known. The best example of 
such a system is a hydrogen atom. Elsewhere we shall 
present the K-R functions for different energy levels of 
the hydrogen atom [nsl . 



C. The Cohen distribution functions 

We have already emphasized the marginal properties 
of a phase space distribution function. An interesting 
question arises about a general form of the distribution 
with the correct marginals. T his p roblem has been posed 
and solved by Cohen in 1966 [|6[ 0. The most general 
distribution P[q,p) with the proper marginals has the 
form of a double Fourier transform of a function 



A{q\p') = J e-^P'i/^ ^*{^ - q'/2) ^{^ + q'/2) 



d^ 



(12) 
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multiplied by an arbitrary function ^{q',p') satisfying 
the relations 



<f((z',o) = i = <i>(o,p')- 



(13) 



In the literature devoted to optical processing of classi- 
cal signals, the function A{q',p') is called the Ambiguity 
function. 

Therefore, the Cohen joint distribution functions la- 
beled by functions $ are given by the following equation 



P^{q,p)=J^[^A] 
1 



(14) 



e'<-P'<J-<j'py''^{q' ,p')A{q' ,p')dp'dq' , 



where by we have denoted a double Fourier transform. 

The Wigner distribution function is obtained by sub- 
stituting ^{q',p') — 1 in the Cohen distribution func- 
tions formula. The K-R distribution function is obtained 
for <i>{q',p') = eM-ip'q'/{2h)]; ^q',p') = cos[p' q' / i2h)] 
leads to the Margenau-Hill distribution. 



III. WIGNER- WEYL TRANSFORMATION 
A. The {q,p) phase space description 

The Wigner- Weyl association of classical phase space 
functions A{q.,p) with quantum operators follows form 
the property: 

{^(g,p)}ordcring = (15) 

dqdpA{q,p) {5{q - q)5{p - p)}ordcrmg- 



The problem of defining a quantum operator in phase 
space lies in the ordering of operators q and p. As an ex- 
ample of such association we shall take a phase space 
density distribution, whose quantum average leads to 
the density operator of the system. Using the Fourier 
decomposition of the Dirac delta functions we see that 
the ordering becomes equivalent to an ordering of the 
Heisenberg-Weyl algebra operators: 



{5{q ~ q)5{p - _p)}ordorm, 

1 



(16) 



dq dp e " ief^e 

I J ordc 



This formula shows that there is no natural unique gen- 
eralization of the classical probability density in phase 
space, because it is possible to have arbitrary classes of 
operator orderings. It has been recognized by Wigner 
and Weyl that various quantum distribution functions 
can be associated with different orderings of operators 
q and p. For example, the Wigner distribution function 
corresponds to Wigner- Weyl ordering, which is obtained 
by putting q and p operators in the same exponent in Eq. 
(p^). We shall show that the K-R distribution function 



corresponds to a special ordering called the anti-standard 
ordering. 

One can easily find that the Ambiguity function, 
Eq. (|l2|) , may be expressed also as 



A{q',p') 



1^) 



(17) 



Using the above equation, the definition of the K-R dis- 
tribution function can be formulated as follows 



K{q,p) 



1 



(27rft)2 



,1,1 i(p'ii-q'p} i 

dp dq e " 



By comparison with Eq. ( |16D we find that the Kirk- 
wood distribution function corresponds to an ordering 
such that all p operators are on the left of all q opera- 
tors. Such an arrangement of the canonical operators is 
called the anti-standard ordering. Analogously one can 
find that the complex conjugation of the K-R function 
corresponds to the standard ordering (all operators q are 
on the left followed by all p operators) . The real part of 
K-R distribution function is associated with a symmetric 
superposition of the anti-standard and the standard or- 
dering, i.e. \{q^p"^ +p"^q^). Note that such an ordering 
is not equivalent to the Wigner- Weyl ordering leading to 
the Wigner function. 



B. Coherent state phase space description 

So far various properties and definitions of phase space 
distribution functions were formulated in position or mo- 
mentum representations. Such a parameterization seems 
to be the most natural to study phase space properties 
of particles. However, Glauber and Cahill have pointed 
out that the coherent state representation is more natu- 
ral and useful while dealing with phase space distribution 
functions describing the quantum states of light [|| . 

We shall now use the following notation: Greek letters 
(a, (3, etc.) designate complex variables; a, annihi- 
lation and creation operators; t>{a) — exp(Q!a^ — a*d) 
denotes displacement operator; and all integrations are 
taken over the whole complex plane. Using the formal- 
ism of coherent states Glauber and Cahill have shown 
that it is possible to define an s-parameterized class of 
distribution functions simply related to the Wigner dis- 
tribution function. These quasi-distributions are defined 
as a complex Fourier transform J- of the s-ordered char- 
acteristic function: 



defined as 
C(/3,s) ^e^l/'I'Tr 



d^(3 



e°''^*^°'*'^C{p,s), (18) 



- /3at -/3*a 



e^"'l Tr[/5D(/3)]. 



(19) 



In the equation above p is the density operator of the 
investigated system. The normally ordered form of D(/3) 
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allows to perform the integral explicitly, and then Eq. 
(\wi) can be rewritten as 



W{a,s) = 



-Tr 



7r(l - s) 

with the operator n(s) defined as 



/5D(a)n(s)Dt(a) (20) 



Il{s) 



s + 1 
s - 1 



(21) 



The continuous parameter s (which is required to be real 
and to satisfy an inequality s < 1) corresponds to differ- 
ing ordering of the creation and annihilation operators. 
Three values: s = 1, s = 0, and s — ~1 correspond to 
normal, symmetric, and anti-normal ordering, that lead 
to the Glauber-Sudarshan P-representation, the Wigner 
function, and the Q-representation, respectively. 

These results have been generalized by Agarwal and 
Wolf in They have proposed the following general 

formula for the quasi-distribution functions: 



(22) 



where fl{f3,P*) is an analytic function of the complex 
variable P that has no zeros. The condition 0(0,0) = 1 
guarantees the normalization 



(23) 



Obviously, such functions do not provide automatically 
the correct quantum marginals, and only for very selected 
functions from the class given by Eq. (^) , we can satisfy 
the relations from Eq. (M). 

In the following Section we shall show that the K- 
R distribution corresponds to a specific form of the 
function. Moreover, with the help of this formalism we 
shall define a new class of a-ordered K-R distribution 
functions in full analogy to the s-ordered functions of 
Cahill and Glauber. 



IV. GENERALIZED K-R FUNCTION AND 
QUADRATIC ORDERING 

A. Definition of generalized K— R distribution 

A new class of quasi-distributions involving a quadratic 
ordering function: 0(/3, /?*) = exp ((t(/9* ^ — /3^)/4) is de- 
fined as a Fourier transform 



C(/?,a), (24) 



of the characteristic function given by 



C(/3,a) =e^^^Tr(pD(/3)), 



(25) 



The continuous parameter a is real and there are no lim- 
its to its value. The K-R distribution function corre- 
sponds to = 1, for (T = we again obtain the Wigner 
distribution function. 

As in the case of s-ordered Cahill-Glauber distribution 
functions, we can perform the integrals and then rewrite 
Eq. (H) as 



K{a,a) = 



ttVT 



=Tr 



where 



K(ct) =: e" 



(26) 



(27) 



In the above formula : : denotes normal ordering of 
operators a, . Using the identity i"^" — : e^*^^)"^" : we 
can transform Eq. (|27|) into 



Kia) 



e 1+° 



1 



e 1+" 



(28) 



Equations (|26[-|27|) define a cr-ordered class of gener- 
alized K-R phase space distributions. From the prop- 
erty K'^{a,a) = K{a,—(T), it is clear that the only real 
function in this class is the Wigner function which corre- 
sponds to cr = 0. 

The operator K((t) is not Hermitian but has a finite 
trace 



Tr 



K(a) 



TT / — 



Let us define operators IC{a,a) as 

^(a,cr) = D(a)K(CT)D^(a), 



(29) 



(30) 



and study its properties. We find that sets IC{a, 1) and 
JC^{a, 1) form a complete, orthogonal basis systems with 
respect to the scalar product 



Tr 



/Ct(a,l)/C(/3,1)J =-<5(a-/3). (31) 



That means that the K-R phase space distribution, up to 
a normalization factor, is an expansion coefficient of den- 
sity operator in the basis /C^ (a, 1). An arbitrary operator 
p may be expanded as 



p^V2 j dPaK{a,l)lC\a,l), 
or, equivalently, 
P 



= \/2 y <faK{a,-l)]C{a,l). 



(32) 



(33) 



From Eq. (32) is clear that knowledge of K-R distribu- 
tion of the state p is equivalent to the knowledge of the 
state itself. The same holds for the complex conjugate of 
K-R distribution function. 
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Connection with squeezed states 



From the formula given by Eq. (g8|) we see that the 
K-R function corresponds to a special case when the op- 



erator 



becomes a projector, as: 



lim , ^ 



(34) 



Thus, as regards the K-R distribution, the formula given 
by Eq. is reduced to 



K{a,l) 



V2 



Tr 



pD(a)e^|0)(0|e-^D^(a) . (35) 



One can easily find that the real part of the K-R dis- 
tribution, up to a normalization factor, corresponds to 
a projection of the density operator on the off-diagonal 
elements of the density matrix of the linear superposi- 
tion given by |$). This connection with the coherent 
squeezed states projection points on an important prop- 
erty of the K-R function, that will be particularly useful 
in the phase space visualization of various properties of 
squeezed states. The relation of the K-R distribution 
to a projection on squeezed states, gives an operational 
meaning of such a phase space distribution. The K-R 
phase space function of a given quantum state at the 
phase space point a is just a projection of this state on 
the off-diagonal elements of the density matrix of the su- 
perposition of infinitely squeezed coherent states. 



The above equation may be rewritten in more intuitive 
form if we refer to the definition of the squeezed state . 
The general coherent squeezed state |a, £,) = E)(q;)S'(^)|0) 
is obtained by a combined action of the displacement and 
the squeezing operator 



(36) 



on the vacuum state. The squeezing operator can be 
written in an ordered form 



expl--a 



t2 



a'a+4 



exp ( -a- 



(37) 



where fi and v are functions of squeezing parameter ^ 
given by: /i = cosh|^| and v = e"^« sinh|f|. Using these 
relations we find that for v / ji = 1 and (j)^ = Q: 



D(a)e'^ |0) = Vcoshr |a, -E), 



(38) 



where ^ = T — arctanh(l) = oo. This means, that left- 
hand-side of the equation above corresponds to the in- 
finitely squeezed state. Hence, the definition of the K-R 
function given by Eq. ( |35| ) takes the form 

j^(«,l)= ^^°'^^ Tr[p|a,-F)(F,a|]. (39) 

TT 

In the similar way we can represent the real part of K-R 
distribution as 

Re[if(a,l)] = ^^^^^il Tr [p(|a, -E) (E, a| 



2-K 



+ |a,E)(-E,a|)]. (40) 



This formula provides a physical interpretation of the 
K-R distribution function in terms of a projection of the 
density operator into a combination of squeezed states. 
Let us consider a linear superposition of two infinitely 
squeezed coherent states 



1$) =AA^(|a,-E> + |a,E)). 



(41) 



V. QUANTUM INTERFERENCE IN PHASE 
SPACE 

All phase space functions from the Cohen class of dis- 
tributions, Eq. (p^), are bilinear in 'J. This provides 
a transparent exhibition of quantum interference. Eor a 
linear superposition of quantum states 



I*) 



l*i> 



1*2 



the corresponding phase space distribution function takes 
the form 

P = Pl+P2+P^nt, 

where Pi, P2 correspond to distribution functions of 
states Y^i) and |^'2), respectively, and Pint denotes the 
interference term. We shall study the properties of the 
interference term for the superposition of two plane waves 
with momenta pi and p2 {h = I): 



~ exp(zpig) + exp{ip2q) 



(42) 



Based on these waves, we shall present a simple argu- 
ment indicating the differences between the quantum in- 
terference patterns exhibited with the help of the Wigner 
function and the K-R function. 

The location of the interference terms of the Wigner 
function is well known. The Wigner function for such a 
superposition is 



Wiq,p) 



5{p~pi) + 5{p 
2 cos {qisp) 5 {p - 



-P2) 



(43) 



where Ap = P2—P1 and p — ^ 



Ei±P2. , The characteristic in- 
terference term cos (qAp) is isolated in phase space from 
the classical points described by the two momenta pi,2, 
and is located between these two incoherent terms at the 
mean momentum p. 

The generalized K-R function for such a superposition 
is given by 

K{q,p)^ S{p-pi) + 5ip-p2) (44) 
+ 5{p-p- o-Ap) Ap + s{p~p + a hp) e-*«^. 
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which leads to the following expression for the real part 
of the generalized K-R distribution 

Re[K{q,p)]^6{p-pi)+5{p~p2) (45) 
+ cos (gAp) 5 {p — p — a hp) + 5 {p ~ p + a tsp) 



We see that the interference term still oscillates like 
cos(gAp), but the locations of these oscillations have 
moved to the points different then in VF((7,p), Eq. (43). 



Analyzing formula (45) we find that parameter a shifts 
the interference term along the momentum axis. For 
CT = 0, which corresponds to the Wigner function, in- 
terference term is located at the mean momentum p. For 
(7 = 1 (case of the K-R distribution) these oscillations 
split and shift to the locations defined hy p = pi and 
p — P2 yt^ the phase space. Increase of a beyond unity 
will move interference term further apart, so that they 
not only no longer overlap each other but appear outside 
the physical location of the state. 

The main difference between the quantum interference 
in the Wigner representation and the K-R representation 
is the location of the interference terms. This difference 
follows from the fact that the K-R function depends lo- 
cally on the phase space properties of the wave function, 
while for the Wigner function this relation is nonlocal. 
This is why the oscillations of the Wigner function occur 
at a position in phase space which is different from the 
classical location. 

It may appear, that the K-R function provides a less 
readable representation of quantum interference, because 
the oscillating terms cannot be isolated from the incoher- 
ent location. However, the advantage of the K-R phase 
space representation is its local relationship with the po- 
sition and momentum wave functions. The Wigner dis- 
tribution function can have oscillations at points where 
the wave function is vanishing. From the definition of the 
K-R function we see that this phase space distribution 
has to vanish when the wave functions "^{q) and ^'(p) 
vanish in q and p. 



VI. LINEAR SUPERPOSITIONS IN THE K-R 
REPRESENTATION 

A. Coherent states and Fock states 

As an example of the structures that can be obtained 
from Eq. (|4|) we consider generalized K-R distribution 
function of a coherent state |ao). By substituting p — 
|ao)(ao| into Eq. ( [2^ ) after simple calculations we get 
the following formula for generalized K-R distribution 
function of coherent state lao): 



Kp{a,a) 



7^^/^ 



(46) 



exp 



2\a — QfoP (T{ot — ao) 




,4 



0,2 



-4 -2 Q 2 4 

FIG. 1: The K-R distribution function for the vacuum state 
|0). 



For (T = it reduces to the well-known formula for the 
Wigner distribution function of a coherent state. For 
other values of a the generalized K-R distributions given 
by Eq. (^^ also consist of a Gaussian bell shape centered 
around a point (Re(Q;o), Im(Q;o)) but modified by a phase 
factor 



e — exp 



1 



(Im(Q;) — Im(Q;o)) (Re(a) — Re(ao)) 



The real part of this distribution has an oscillating term 
cos (5, making the quasi-distribution non-positive. For 
(7^0, this oscillating term corresponds to the plane 
wave from Eq. (|^) with the phase <p(g,p) proportional 
to a factor , . With increase of the value of cr, these 
oscillations quickly become more and more rapid, achieve 
maximum frequency for cr = 1 and then slowly vanish. 

Substituting ao = and a = I into the Eq. ( p6| ) we 
obtain formula that describes K-R distribution of the 
vacuum state |0): 



K\o){o\ia, 1) = — exp 



a 
Y 



(47) 



Figure 3 shows the real part of the expression given by 
Eq. (^). In accord with the previous description, it is 
a Gaussian function in position and momentum, modu- 
lated by the plane wave phase factor cos[2Re(Q;)Im(Q;)]. 

Formula for K-R distribution of the one photon Fock 
state |1) is given by 



K\i)ii\{a, 1) 



V2(a2 - a* 



IT 

,2 /^,*\2 



exp [ - y 
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,4 



0,2 



-4 -2 Q 2 4 

FIG. 2: The K-R distribution function for the one photon 
Fock state. 



The real part of this equation is shown in the Fig. ||. 
The difference between expressions for the K-R function 
of the vacuum state and the one photon state is an ad- 
ditional multiplicative amplitude iqp of the oscillating 
phase factor. This amplitude comes from the product of 
the one photon wave functions in momentum and posi- 
tion representations, respectively. 



B. Superposition of coherent states 

Another simple but interesting and instructive exam- 
ple that we shall study is a superposition of two coherent 
states N{\ao) + \ - ao)). For this state, from Eq. (H) 
we derive: 



This is illustrated in Fig. where we have plotted the 
real part of Eq. (|4|) for ao — 3 and several values of 
parameter a. We present the real part of distributions 
rather then the imaginary part (that does not contribute 
to marginals), or the absolute value (in which no oscilla- 
tions appear). In Fig. ^ we present the same functions 
using contour plots. In this representation it is apparent 
how oscillating interference term moves with the change 
of the value of the parameter a. As was the case of 
the superposition of plane waves, for cr = interference 
term is located exactly between the terms that do not 
correspond to the interference. With the increase of the 
value of the parameter a interference terms split away 
and move along Re(Q!) axis, as they are centered around 
phase space points (Re(a), Im(a)) = (±(Tao,0). Again, 
for a — 1 oscillating interference terms overlap the terms 
corresponding to the two individual coherent states. As 
we have already mentioned, the K-R distribution func- 
tion is a special case of a distribution which is always 
equal to zero when wave function is zero at certain point. 

Figure ^ shows the same distribution functions when 
the initial state is rotated in phase space. As an example 
we have chosen previous initial state rotated by <^ = ^ 
in phase space, i.e. ao = 3 . The Wigner function of 
a "rotated state" is simply the rotated Wigner function. 
This property is the basis for tomography of the Wigner 
function. By contrast, for other distributions (with pa- 
rameter (T 7^ 0) this simple relation does not hold, and 
tomography would not work. 

Let us consider behavior of the terms corresponding to 
the individual coherent states in Fig. ^. Although the 
center of such a term is located exactly in this point of 
phase space where coherent state is centered, closer ex- 
amination tell us that whole term is not "rotated" but 
rather "shifted" by appropriate values parallel to the co- 
ordinate axes. It is clearly seen that for cr 7^ gen- 
eralized K-R distribution functions single out p and q 
axes as compared to the other directions in phase space. 
Interference terms are just rotated but in the opposite 
direction. 



2N 



[<t(q 



aS)^-a(a-co)''-2|c~Qo| 

TTr2 



-t- e 



[^(o*+c.g)^-a(a + e.o)^-2|c + Qop] 
a(o*+c*)3-a(o-Q0)2-2(o*+c.g)(c-e»0)] 

[.T(a*-Qg)^-,T(Q + ao)^-2(a*-Qg)(Q + QO)l 



-2|aol' 



,(48) 



where 



N 



2 + 2e 



-2|ao|^ 



Besides terms corresponding to two individual coher- 
ent states there is an interference cross term. Its local- 
ization changes with the change of parameter a in the 
same way as for the superposition of two plane waves. 



C. Superposition of coherent squeezed states 

In the previous section we have connected the real part 
of the K-R function with the expectation value of off- 
diagonal elements of density matrix corresponding to su- 
perposition of two infinitely squeezed coherent states. It 
is interesting to see how K-R distribution evolves if we 
change the squeezing parameter. 

In Fig. II we have plotted the real part of the K-R 
distribution function of superposition of squeezed vac- 
uum states A/'dO + I — 0) for several values of squeezing 
parameter ^. The K-R distribution apparently shows 
that an investigated state consist of two perpendicularly 
squeezed states. Obviously, the closer parameter ^ gets 
to unity, the more significant are effects of squeezing. 
Moreover, also basic properties of K-R distribution are 
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FIG. 3: Generalized K-R distribution function of superpo- 
sition of two coherent states NQao) + \ — ao)). Here we have 
chosen ao = 3, and a) a = (the Wigner distribution func- 
tion), h) a = 0.5, c) a = 1 (the K-R distribution function). 



FIG. 4: Contour plots of generalized K-R distribution func- 
tion of superposition of two coherent states A'^dao) + | — ceo))- 
liere we have chosen ao = 3, and a) a = (the Wigner distri- 
bution function), h) a = 0.5, c) cr = 1 (the K-R distribution 
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C = 0.1 
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FIG. 5: Contour plots of generalized K-R distribution func- 
tion of superposition of two coherent states N{\an) + [ — ao})- 
Here we have chosen qq = 3 -|- i3, and a) a = (the Wigner 
distribution function), b) cr = 0.5, c) cr = 1 (the K-R distri- 
bution function). 




i = 0.5 




^ = 0.7 




^ = 0.9 




FIG. 6: The K R distribution function of superposition of 
two squeezed coherent states Af{\^) + \ — §))• 
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emphasized with the increase of the value of the param- 
eter ^: the oscihating cos[2Re(Q;)Im(a)]-hke structure is 
more and more distinct. As we have mentioned before, 
the K-R distribution single out p a, q axes from other 
phase space directions. Looking at the changes caused 
by the increase of squeezing, one gets easily convinced of 
a close connection between the K-R distribution and a 
superposition of infinitely squeezed coherent states. 

VII. SUMMARY 

We have presented a new class of phase space quasi- 
distribution functions with correct momentum and posi- 
tion marginal properties, that contains the Wigner dis- 
tribution function and the K R distribution as special 
cases. We have shown how quantum interference appears 
in phase space if such functions are used for its investi- 
gation. In particular, we have focused on similarities and 
differences between the Wigner function and the K-R 



distribution. 

We have emphasized the most important properties of 
the K-R distribution function: the fact that this func- 
tion fully characterizes the quantum state; that the K-R 
function corresponds to the anti-standard ordering of q 
and p operators; and that the real part of K-R distri- 
bution is in natural way connected to coherent squeezed 
states. 

Most of the properties of the generalized K R distri- 
bution function we have presented using the formalism 
of coherent states with the full analogy to the Glauber 
and Cahill s-ordered quasi-distributions. 
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